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We investigate a Kondo lattice model with correlated conduction electrons. Within dynamical 
mean-field theory the model maps onto an impurity model where the host has to be determined 
self-consistently. This impurity model can be derived from an Anderson-Hubbard model both by 
equating the low-energy excitations of the impurity and by a canonical transformation. On the level 
of dynamical mean-field theory this establishes the connection of the two lattice models. The impu- 
rity model is studied numerically by an extension of the non-crossing approximation to a two-orbital 
impurity. We find that with decreasing temperature the conduction electrons first form quasipar- 
ticles unaffected by the presence of the lattice of localized spins. Then, reducing the temperature 
further, the particle-hole symmetric model turns into an insulator. The quasiparticle peak in the 
one-particle spectral density splits and a gap opens. The size of the gap increases when the corre- 
lations of the conduction electrons become stronger. These findings are similar to the behavior of 
the Anderson-Hubbard model within dynamical mean-field theory and are obtained with much less 
numerical effort. 
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I. INTRODUCTION 

The usual explanation for the formation of heavy 
fermions in compounds with rare-ejarth or actinide ele- 
ments is based on the Kondo effect .Idfl Thereby, the char- 
acteristic low-energy scale arises from the spin-screening 
of the local moments by a non-interacting electron gas. 
The periodic Anderson model, as well as the Kondo lat- 
tice model, are considered as the most promising candi- 
dates to at least qualitatively describe the rich physics of 
these materials. However, these models faiE to describe 
correctly the temperature scale of the heavy-fermion be- 
havior in the electron-doeed cuprate Nd2-xCexCu04 dis- 
covered a few years ago.Q It has therefore been suggested 
that the strong correlations among the doped carriers 
which are usually neglected in the Kondo effect play a 
crucial role jand may enhance the Kondo temperature 
significantly.^ Subsequently, efforts have been undertaken 
in order to study the influence of the correlations among 
the conduction elect rons_pp_. the Kondo effect, both for 
impurityOlj and latticetjlij models. 

These correlations are usually introduced by adding 
a Hubbard-type interaction among the conduction elec- 
trons to the Anderson model. In Ref. 15 the resulting (pe- 
riodic) Anderson-Hubbard mode.1 wfls investigated within 
the dynamical mean-field theory.E3llII It was shown that it 
maps onto an impurity model where the impurity consists 
of a unit cell of the lattice model, i.e., of two correlated 
orbitals. This model still includes charge fluctuations on 
the / orbital which should not influence the low-energy 
behavior in the Kondo limit. It would be desirable to 
eliminate these degrees of freedom as this reduces the 
number of states of the impurity and thus reduces the 



computational effort to investigate the impurity model. 
Hence, we are aiming at a Kondo-Hubbard model as the 
counterpart of the Anderson-Hubbard modeLuist, as in 
the case of uncorrelated conduction electrons £3 Ej 

In the next section, we introduce the Kondo-Hubbard 
model and map it to an impurity model using the dy- 
namical mean-field approximation. In Sec. [II we relate 



the derived effective impurity model with the one of the 
Anderson-Hubbard pandel. We have applied the resolvent 
perturbation theoryEil to investigate the effective impu- 
rity model for the Kondo-Hubbard model. We present 
some results in Sec. |^ and finally conclude in Sec. 0. 



II. THE KONDO-HUBBARD MODEL 

The (periodic) Kondo-Hubbard model is given by 
Hk = H, + Hs, (1) 

ka ^ ^ i 
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(3) 



Here, the c operators refer to the conduction electrons 
which are described by a Hubbard model with on-site re- 
pulsion Uc- S!j^ denotes the spin of an electron at site i 
in the conduction band orbital c and the impurity orbital 
/, respectively. For simplicity we consider in this paper 
the particle-hole symmetric model only, which implies 
/i = 0. 



We closely follow the treatment of the usual Kondo 
lattice model in infinite dimensions given in Ref. ^. We 
introduce Fermion operators dia to represent the / spin 
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with local constraints 
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for all i. To enforce the constraints we add 
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to the Hamiltonian and take the limit [/ ^ oo.E3 The 
Haniiltonian i/x includes three two-particle parts pro- 
portional to J, Uc and [/, respectively. Given the non- 
interacting Green's function, G'o(fc,Ci;), of the Kondo- 
Hubbard model 
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the self-energies are defined by Dyson's equation 

G{k,u;) = {Go\k,u;) -J:{k,u;)y\ 

The dynamical mean-field theo« assumes a momentum- 
independent self-energy matrixEj I](fc, (jo)-,— > S(ti') which 
holds in the limit of infintc dimensions .cJ Following the 
arguments given in Ref. |l5| we can calculate S(a;) from 
the impurity model 
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Equation (h2|) is valid for a semi-elliptic density of states 
of the conduction electrons 

Poi^) ^ ^ VW^ ~ e^ . (13) 

1 as unit of energy. 
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with the self-consistency equation 
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In all our calculations we choose W 

In order to solve numerically the impurity model (|9|) 
we make use of the resolvent perturbation theory intro- 
duced in Ref. ^. It serves as a basis for the nonjja-ossing 
approximation to the Anderson impurity modeled which 
has been successfully extcijwicd to investigate the finite- [/ 
Anderson impurity model, ^ the Emery model in the dy- 
namical mean field theoryJ£3 aad the Anderson-Hubbard 
model in this approximation.ll3 As described in Ref. UB, 
the occuring coupled integral equations are-solved numer- 



ically by introducing defect propagatpsstSi and making 



use of the fast Fourier transformation.! 



III. RELATION WITH THE 
ANDERSON-HUBBARD MODEL 

It is well known that the Anderson and Kondo impu- 
rity mojdei-are mutually related,ll2 as well as the lattice 
models.Ej€j We expect that also the lattice models which 
include interactions among the conduction electrons show 
the same low-energy behavior. Within the dynamical 
mean-field theory this should be true provided the ex- 
citation energies of the corresponding impurities arc the 
same in both models. In Ref. Ha we have derived the 



effective impurity model Hp^i 



H.r 



Hioc for the 



Anderson-Hubbard model. The local part reads 



^loc = -^f^c E 4c. + U^n'^nl + VY, Hfa + H.c.) 
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Note that the mapping being discussed here is based on 
the impurity models and not on the original lattice mod- 
els. Therefore, the same restrictions apply as to the dy- 
namical mean-field theory or the limit of infinite coordi- 
nation number of the underlying lattice. In particular, 
spatial correlations are not, included, so that an RKKY 
interaction will not occur p3 in contrast to the treatment 
in finite dimensions. E2I 

It has been found in Ref. ^ that two characteristic 
features occur at low energies for the Anderson-Hubbard 
model. At intermediate temperatures, a quasiparticle 
peak is formed at the chemical potential which corre- 
sponds to thc-Quasiparticle peak observed for a pure Hub- 
bard model n3'BE2l In this temperature regime the pres- 
ence of the / orbitals does not influence the low-energy 
behavior. When lowering the temperature a gap opens. 
As in the case of free conduction electrons, this is inter- 
preted as a hybridization gap which arises from the level 
crossing of the quasiparticle band with dynamically cre- 
ated states-, at the chemical potential (Kondo-resonance 
states) .E3E3 

Both features are related to transition energies of the 
impurity which generate a Kondo effect on a certain en- 
ergy scale when coupled to a bath. The eigenenergies 
of the impurity ( |l4| ) where the / orbital is singly occu- 
pied when V^ — > are listed in Tab. |[ The Hubbard-like 
quasiparticle peak is not influenced by the presence of 
the / orbitals and results therefore from the c charge 
degree of freedom. The relevant transition in ifioc is 
/^c" — > f^c^{S) {S stands for singlet) with its energy 
given by AEi = E{fc°) - E{fc\S)). The opening 
of the gap is found at lower temperatures and it is re- 
lated to the singlet-triplet excitation f^c^{T) -^ f^c^{S) 
(the lowest transition energy of the impurity) with energy 
AE2^Eifc\T))-Eifc\S)). 

The two impurity models generate the same low-energy 
behavior, provided their corresponding low-energy tran- 
sitions have identical energies. The eigenenergies of 
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TABLE I. Eigenenergies of the impurity, Hioc, in the An- 
derson- (AHM) and Kondo-Hubbard model (KHM) in the 
symmetric case. The indicated occupation is the one for 
V ^ 0. S{T) stands for singlett (triplett) configuration, resp. 
States not included for the AHM have zero / occupation in 
the limit V —^ and no counterpart in the KHM. 



the Kondo impurity (|10|) are listed in Tab. | as well. 
We determine the parameters of the Kondo-Hubbard 
model by equating the corresponding transition ener- 
gies. The exchange coupling constant J is obtained from 

Eifc^iT)) - E{pc^{S)) ^ AE2 ^ J 
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The Hubbard interaction in the Kondo-Hubbard model, 
which we denote by Uc in order to distinguish it from the 
Anderson-Hubbard model, is obtained from E{f^c^) — 
E{fc^S)) = AEi = Uc/2 + 3/4 J 



Uc = g(C// 




(16) 



Whereas the appearance of the exchange coupling J as 
well as its magnitude for small V is expected (the inter- 
mediate state requires a doubly occupied c orbital with 
energy Uc/2), the modification of the local interaction of 
the conduction electrons might be surprising. However, 
if / charge fluctuations are possible, a doubly occupied c 
orbital may lower its energy by virtual transitions onto 
the / orbital. To compensate, the effective repulsion is 
lowered when the / occupancy is fixed. 

The same relations between the parameters of the 
Anderson- and Kondo-Hubbard model are obtained by 
performing a canonical transformation of the molecule 
where one eliminates the charge degrees of freedom of 
the / orbital in perturbation theory in the hybridization 
V. This is shown in App. |^. 



T = 0.025 




FIG. 1. Local density of states Ac(z) in the 
Kondo-Hubbard model for Uc = 1 and J — (solid), 0.05 
(dotted), 0.1 (dashed), and 0.15 (dot-dashed line) at different 
temperatures. 



IV. RESULTS 

In Fig. Q we display the local spectral functions of the 
conduction electrons which were obtained by the treat- 
ment outlined in Sec. y for Uc ^ 1 and various J at 
different temperatures T. The results found are typical 
for the Kondo-Hubbard model. At comparatively high 
temperatures (T > 0.5), only two broad maxima are seen 
which correspond to upper and lower Hubbard band of 
the conduction electrons. When lowering the temper- 
ature, a peak emerges at the chemical potential. The 
spectra at these temperatures are almost independent of 
the coupling J to the / spin for small J < 0.1, and co- 
incide with the spectra calculated for the pure Hubbard 
model (J = 0). Therefore, we attribute thia,pfiak to the 
quasiparticle peak of the Hubbard modelE300 which 
describes the conduction electrons. At these moderate 
temperatures the c and / electrons are almost decou- 
pled. The quasiparticle peak corresponds to the c^ — > c" 
transition in the c-f molecule. Coupled to the bath, this 
transition gives rise to a Kondo effect, which sets a char- 
acteristic energy scale. Assuming J = and a structure- 
less bath with A(0) = 2/7r, the impurity model reduces 
to the symmetric Anderson model jand the Kondo tem- 
perature is estimated to Tk '^ 0.5E3 (This rather large 
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FIG. 2. Gaps of the Kondo-Hubbard model for U^ = 1.0 
and 1.5. The hnes are guides to the eye. 



scale is due to the large coupling of the impurity to the 
bath which is of the order of the band-width D.) No- 
tice that the given value of Tk provides only a rough 
estimation of the corresponding energy scale in our cal- 
culation. First of all, the assumption of a structureless A 
is not valid. Moreover, neglecting vertex corrections in 
the NCA for the symmetric Anderson model, aa^we have 
done, underestimates the Kondo temperature.o 

When we reduce the temperature further, we observe 
that the quasiparticle peak gets split (Figs. ||). This 
splitting results from the influence of the / spins since 
its characteristic temperature, as well as its magnitude 
depend on J. Therefore, the / and c subsystems are no 
longer decoupled in this temperature regime. At the low- 
est temperatures that we reach numerically, and which 
are lower when J becomes smaller, a gap begins to open 
as is expected for the symmetric model. Qualitatively, 
these findings fit to the scenario encountered in the An- 
derson model with free conduction electrons, where below 
a characteristic temperature a local Kondo effect takes 
place, in which the conduction electrons screen the / 
spin, which leads to a resonance at the chemical poten- 
tial in the particle-hole symmetric case. These dynami- 
cally generated states cross the conduction-band states, 
one finds a splitting of the conduction band with a gap 
at the ch emi cal, potential, and the system becomes an 
insulator.k3Eay 

Quantitatively, however, these effects depend not only 
on the exchange interaction J but also on the local Hub- 
bard interaction of the conduction electrons Uc- This is 
seen in Fig. where we plot the gaps measured as the 
distance between the two maxima for C/c = 1 and LS. 
vs. J. As discussed for the Anderson-Hubbard model,ta 
our numerical procedure does not allow to proceed to 
very low temperatures. As the gaps have not yet com- 
pletely converged we only repeat the following qualita- 
tive remark which applied also for the Anderson-Hubbard 
model: The larger Uc, the larger the resulting gap, and as 



in the case of the Anderson-Hubbard model, the quasi- 
particles of the Hubbard model are responsible for the 
screening of the / moment. 

Next, we compare thesci-£ndings to the results of the 
Anderson-Hubbard modcO in Fig. (left). The pa- 
rameters of the Kondo-Hubbard model, Uc = 1.479 and 
J = 0.04851, were chosen such that the model origi- 
nates from an Anderson-Hubbard model with Uc = 1.5, 



5, and y = 0.2 — see Sec. HI 



The most obvious 
difference is the presence of the / peak at z ~ zLUf/2 in 
the latter case which corresponds to the possible charge 
fluctuation on the / orbital. Except for this, we find that 
the two spectra are similar. Since the Kondo-Hubbard 
model was derived from the Anderson-Hubbard model by 
projecting out charge fiuctuations on the / orbital this 
should be expected for a small hybridization V where 
these charge fluctuations are unimportant. 

However, close to the chemical potential the two spec- 
tra only qualitatively agree. This should not be too 
surprising, as in the canonical transformation the one- 
particle operators should be tcansformed as well in order 
to obtain equivalent spectra,E3 see also App. ^. Besides 
higher Green's functions, one finds contributions from 
Qcf and Gf of the Anderson-Hubbard model to the c 
spectral function of the Kondo-Hubbard model. We end 
this section by reporting a strange observation. If we ig- 
nore the change of Uc when passing from the Anderson- 
Hubbard model to the Kondo-Hubbard model, a differ- 
ence between the two spectra can hardly be seen at small 
V, Fig. ^ (right). The rather strong influence of Uc on the 
spectra (C/c changed by 1.4% from Fig. || (left) to Fig. ^ 
(right)) is due to the proximity to the metal-insulator of 
the bare conduction electrons which in our scheme occurs 
at [/c ~ 1-7. For larger values of V this finding does not 
hold true, but the agreement in total gets worse since 
charge fluctuations become more important. 



V. CONCLUSIONS 

In summary, we investigated a Kondo lattice model 
with correlated conduction electrons within the dynami- 
cal mean-field theory. To this end we mapped the model 
onto an impurity model in which a unit cell of the lat- 
tice is embedded self-consistently into a bath of free elec- 
trons. This impurity model was studied numerically 
within the resolvent perturbation theory by applying the 
non-crossing approximation. 

The parameters of the Kondo-Hubbard model were re- 
lated to the parameters of the Anderson-Hubbard model 
by matching the transition energies within a single unit 
cell. It was shown that this is equivalent to performing 
a canonical transformation within a single unit cell. 

The resulting spectral functions (local density of 
states) behave similarly as in the Anderson-Hubbard 
model: At relatively high temperatures the conduction 
electrons are almost decoupled from the / electrons and. 
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FIG. 3. Local density of states Ac{z) in the Kondo-Hubbard model (solid lines, J = 0.04851, Uc ~ 1.479, left, and Uc ~ 1.5, 
right) and the Anderson-Hubbard model (dashed lines, Uc = 1.5, Uf = 5,V = 0.2) at different temperatures. 



as in the pure Hubbard model, we observe a quasipar- 
ticle peak emerging when the temperature is lowered. 
When the temperatures is reduced further, the quasipar- 
ticle peak splits and a gap opens. As in the case of free 
conduction electrons, one might envisage that a (local) 
Kondo effect generates resonances at the chemical poten- 
tial which hybridize with the quasiparticle states leading 
to a gap. 

The correlations among the conduction electrons, how- 
ever, do not only show up in the formation of a quasi- 
particle peak at high temperatures, but also influence 
the size of the gap at low temperatures. In agreement 
with the Anderson-Hubbard model, we find that the 
low-energy scale is increased as the correlations become 
stronger. 

Both features result from certain transitions in the im- 
purity model. The first one, which corresponds to the 
Hubbard quasiparticle peak, is related to the (local) tran- 
sition where the c orbital changes from single to dou- 
ble (empty) occupancy, whereas the second one, leading 
to the gap, stems from the singlet-triplet transition of 
the c-f molecule which has a lower energy. The consid- 
ered impurity model is a "minimal" model which contains 
two different excitation energies and the Kondo-Hubbard 
model is a minimal model to generate the two features 
mentioned above. Therefore it deserves continued inter- 
est. 



APPENDIX A: CANONICAL 
TRANSFORMATION 

The traditional way of deriving Kondo models from 
Anderson models is to perform a canonical transforma- 
tion 



Hk = e^i^Ae^ 



(Al) 



such that there is no c-jOiybridization in iJx and the / 
occupancy is conservedJlj Usually this is done in pertur- 
bation theory with respect to the c-f hybridization V. 
We will perform this transformation on the level of the 
impurity models and rewrite our impurity model as 



^A,iinp = Ho + Hi , 

Ho = Ha + Hac + Hc + Hf , 

Hi = Hcf ■ 
The transformed Hamiltonian is given by 
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where SA = [S', A]_ and 5^ ^ V^. Charge fluctuations 
on the / orbital are ehminated by choosing S such that 
there are no contributions of odd power in V. Hence, 

SiHo = -Hi ^ Si = -^Hi , (A7) 

S3H0 = -^SfHi ^83 = —S^Hi . (A8) 

where the Liouvillean Cq is defined by CqA ~ [Ho,A]^. 
The efi^ective Hamiltonian which conserves the / occupa- 
tion is then given by 



i?K,imp = Ho + -SiHi + -SfHi 
Z o 



(A9) 



In order to reproduce the previous results from Sec . Ill 
we will neglect Hac and thus replace Cq in Eqs. ( [A^ ) by 
Cc + Cf. The Liouvillean is easily inverted by decompos- 



ing c^ 



^ -|- cj. where 



cinl 



(AlO) 



are eigenoperators of Cc with eigenvalues —Uc/2 and 
C/c/2, resp. (similarly for /J). We immediately obtain 
^i and find as effective Hamiltonian within the space 
n/ = 1: 



-ffK.imp — Ha + Hac + H^ 
+ JSfSc + SUc 

where we dropped a constant and 
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Uf -Uc Uf + Uc 
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(Uf - U.f {Uf + Ucf 



(All) 
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(A13) 
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These are just the results ( p^ and (EBI) found in the 
previous section by comparing the excitation energies of 
the impurity in both models. 

Note that we neglected the coupling of the c-f molecule 
to the bath in this derivation. In reality the a oper- 
ators should change in the canonical transformation as 
well. This approximation corresponds to just nratching 
the molecular excitation energies in the previous sec- 
tion. Note also that the canonical transformation was 
performed on the impurity model. Therefore no RKKY 
interaction is generated to forth order in V . 

When comparing spectral functions of models related 
by canonical transformations one should bear in mind 
that|-the creation operators should be transformed as 
well.Ea This is seen as follows: The Hamiltonians in both 



models are related by ( |A5[) , the ground states via 



\^K) = e'^\^A) 



(A15) 



Thus, the c correlation function, e.g., of the Kondo- 
Hubbard model is given by 



Ck{z) 






1 



^"z-Hk 



* 



K 



(e^)tc.(e-^) 



S\'i. 



z-Ha 



e-^cle^ 



*a) (A16) 



(and correspondingly the Green's function). In particu- 
lar, the c Green's function in the Kondo-Hubbard model 
is not identical to the one in the Anderson-Hubbard 
model. Given 5*1 wc find to first order 



e-^cie^ ^ci + 2V ■ 



2Uc 



Uf - ui 



-fl\ ■ (A17) 
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